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INTRODUCTION 


The  Astronautics  and  Geodesy  Division  in  conjunction  with  the  Advanced  Projects  Division  of 
NSWC,  with  sponsorship  by  the  Defense  Mapping  Agency,  have  developed  a  receiving  system  that 
allows  the  NAVSTAR  Global  Positioning  System  (GPS)  to  be  used  for  geodetic  point  positioning. 
Previous  experience  with  the  Navy  Navigation  Satellite  System  has  shown  that  observations  of  Dop¬ 
pler  frequency  from  many  satellite  passes  can  be  used  to  obtain  positions  at  the  receiving  site  with  an 
error  of  only  a  few  meters.  It  was  hoped  that  similar  techniques  applied  to  the  GPS  constellation 
would  allow  even  better  results  to  be  realized.  It  was  to  demonstrate  this  potential  that  the  NGRS 
was  developed. 


SYSTEM  DEFINITION 


The  heart  of  the  NGRS  is  the  GPS  receiver,  which  was  developed  by  the  Stanford  Telecom¬ 
munications,  Incorporated,  under  contract  to  NSWC.  The  receiver  was  delivered  in  September 
1978.  Integration  of  the  receiver  and  the  other  equipment  that  comprise  the  NGRS  progressed 
through  the  end  of  1978.  First  data  were  processed  in  January  1979,  and  the  system  was  declared 
operational  at  the  end  of  that  month.  The  truck  containing  the  NGRS  was  driven  to  Yuma  Proving 
Ground  in  early  February  for  the  initial  tests  alongside  the  Mobile  Test  Van,  which  contained  a 
Magnavox  “X-set”  configured  for  static  real-time  positioning.  The  NGRS  remained  at  the  Yuma 
site  until  day  80  of  1979. 


Figure  1  shows  a  block  diagram  of  the  NGRS,  and  a  list  of  the  equipment  is  given  in  Table  1 .  A 
description  of  the  satellites  and  the  signal  structure  is  described  in  a  special  issue  of  the  Journal  of 
the  Institute  of  Navigation1 .  Each  satellite  transmits  on  two  L  band  frequencies:  L,  =  1575.42  MHz 
and  Li  =  1227.6  MHz.  The  bandwidth  of  the  signals  is  about  20  MHz  and  consists  of  a  carrier 
biphase  modulated  by  two  pseudo-random  sequences  having  frequencies  of  1.023  MHz  {C/A)  and 
10.23  MHz  ( P ).  A  data  bit  stream  of  frequency  50  Hz  is  added  to  each  of  the  pseudo-random  se¬ 
quences.  A  replica  of  the  C/A  and  Pcode  is  generated  in  the  receiver  and  correlated  with  the  incom¬ 
ing  signal.  The  replica  is  shifted  in  time  until  a  correlation  peak  is  obtained.  If  this  correlation  peak 
is  maintained,  the  local  code  is  then  locked  to  the  incoming  code  from  the  satellite.  The  dependence 
of  this  correlation  peak  on  bandwidth  is  described  in  Appendix  A. 


The  pseudo-range  measurement  is  obtained  by  recording  the  time  interval  between  the  local 
time  epoch  and  the  same  epoch  received  from  the  satellite.  This  measurement  consists  of  two  parts: 
the  propagation  delay  and  the  time  offset  between  the  local  and  satellite  clocks.  Figure  2  shows  this 
relationship. 


The  analog  Doppler  signal  derived  from  the  reconstructed  carrier  in  the  receiver  is  used  to  ob¬ 
tain  Doppler  counts  on  both  L,  and  L1.  The  time  interval  over  which  the  count  is  made  is  deter¬ 
mined  by  the  local  clock  and  is  nominally  60-s  duration.  A  continuous  running  counter  is  strobed  at 
the  end  of  the  interval  to  provide  the  integer  count.  The  time  from  the  end  of  the  interval  until  the 
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ro  =  EPOCH  OFFSET  (SATELLITE  -  LOCAL) 

r„  -  OBSERVATION  OF  PSEUDO*RANGE  =  t0  +  rr 
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A  VALUE  FOR  r,  IS  FOUND  IN  THE  SATELLITE  DATA  MESSAGE 


Figure  2.  Relationship  Between  the  Time 
Systems  and  Pseudo  Range 


LOCAL  TIME  tM  tj 
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ZERO  DOPPLER  FREQUENCY  =  28750  Hi 


Figure  3.  Measurements  That  Comprise  the  Doppler  Observations 


next  positive  zero  crossing  of  the  analog  Doppler  signal  is  also  recorded.  This  allows  a  precise 
measurement  of  the  true  time  interval  to  be  made.  A  pictorial  description  of  the  process  is  shown  in 
Figure  3. 


The  system  operation  is  controlled  by  a  microprocessor.  To  begin  observing,  the  operator  is 
asked  to  key  in  the  options  he  wishes  to  use  for  the  current  session.  Each  response  is  displayed  on  the 
video  terminal  and  sets  flags  in  the  controlling  program.  After  a  cassette  tape  is  loaded  into  the  tape 
drive,  the  system  is  ready  to  begin  observations.  A  search  is  initiated  for  the  first  satellite  using  the 
Doppler  estimate,  which  has  been  entered  by  the  operator. 


Two  methods  of  operation  are  available:  local  and  remote.  In  the  local  mode,  the  operator 
loads  the  acquisition  information  directly  into  the  receiver  via  the  front  panel  switches.  This  infor¬ 
mation  consists  of  the  satellite  code  identification  number,  the  expected  Doppler  frequency  with  a 
tolerance  of  ±500  Hz,  and  the  differential  ionospheric  delay  between  the  L ,  and  L2  channels.  Ac¬ 
quisition  of  different  satellites  is  made  by  entering  new  information  when  desired.  The  remote  mode 
offers  the  option  of  entering  a  list  of  acquisition  information  and  start  times  into  the  microprocessor 
via  the  the  keyboard.  The  computer  will  then  operate  the  receiver  remotely,  sending  the  new  acquisi¬ 
tion  information  at  the  preset  times. 


When  the  receiver  has  acquired  the  first  satellite,  pseudo-range  measurements  obtained  from' 
the  time-interval  counter  are  sent  to  the  microprocessor  each  6s.  At  the  end  of  a  minute,  as  deter¬ 
mined  by  the  local  clock,  the  two  Doppler  counts  and  the  two  fractional  counts  are  sent  to  the 
microprocessor.  All  the  data  accumulated  during  that  minute  are  then  recorded  on  the  tape 
recorder.  The  system  continues  tracking  in  this  1 -minute  cycle  until  the  satellite  goes  out  of  view,  or 
the  time  arrives  to  search  for  a  different  satellite.  The  tape  holds  about  12  hr  of  continuous  observa¬ 
tions,  and  it  can  switch  to  a  second  drive  when  the  first  is  full.  The  data  from  these  cassettes  are 
transferred  to  a  nine-track  tape  and  are  then  preprocessed.  The  preprocessed  data  are  kept  on  disk 
packs  for  further  processing  as  desired. 


TIME  CORRECTIONS 


Time  corrections  are  necessary  so  that  the  Doppler  observations  and  the  satellite  reference  tra¬ 
jectory  are  referred  to  the  same  time  system.  For  our  purposes,  time  is  kept  at  the  Master  Control 
Station  (MCS)  at  Vandenburg  Air  Force  Base.  This  is  GPS  time  and  is  allowed  to  follow  its  own 
course;  that  is,  it  is  not  corrected  for  the  variations  in  the  earth’s  rotation  as  is  Universal  Time. 


Each  GPS  satellite  has  an  independent  frequency  standard  that  must  be  corrected  to  GPS  time. 
The  corrections  for  each  satellite  are  predicted  by  the  MCS  and  uploaded.  Each  satellite,  in  turn, 
transmits  these  predicted  corrections  to  its  apparent  time  in  its  navigation  message.  The  user  must 
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then  apply  these  corrections  when  he  uses  the  data.  In  this  way,  all  the  data  can  be  converted  to  a 
common  time  system. 


In  the  NCiRS,  the  single  frequency  (L,)  pseudo-range  measurement  is  used  to  relate  the  local 
time  to  GPS  time  (Appendix  B).  The  pseudo-range  observation  gets  its  name  because  it  is  just  a 
t ime-of -flight  range  measurement  treated  as  if  the  clocks  on  both  ends  were  synchronized.  In  reality 
they  are  not  synchronized;  and  so  the  measurement  is  not  a  true  range,  but  includes  the  clock  errors 
plus  propagation  effects.  Figure  2  illustrates  the  measurement.  Epoch  is  an  arbitrary  instant  at 
which  each  clock  registers  the  same  numerical  time.  The  observable  is  T„and  is  time  tagged  by  the 
satellite  epoch.  However,  the  satellite  epoch  is  offset  from  GPS  time  by  t,  (note  that  t,,  as  well  as  t, 
and  t„,  are  not  constant  but  vary  with  time).  The  navigation  message  from  each  GPS  satellite  gives 
the  predicted  values  for  t,  as  a  function  of  GPS  time.  Thus,  the  true  time  (GPS  time)  tag  for  the 
observation  is  the  satellite  epoch  time  minus  t,.  This  is  the  correction  that  needs  to  be  made  to  the 
pseudo-range  time  tags  to  indicate  the  true  time  of  transmission.  With  the  use  of  the  t„  observations, 
the  local  clock  offset  r0  can  be  obtained  by  subtracting  out  the  propagation  delay  t,.  To  calculate  t,  , 
the  satellite  position  and  the  receiver  position  are  needed.  However,  extreme  position  accuracy  is  not 
necess”- '  because  the  relative  velocity  is  much  less  than  the  propagation  velocity. 


If  a  worst -case  relative  velocity  of  1000  m/s  is  assumed,  the  slant  range  changes  by  only  0.001  m 
in  1  ^s.  Light  travels  300  m  in  1  /^s.  Therefore,  a  position  error  of  300  m  in  the  slant  range  would  pro¬ 
duce  only  1-ps  error  in  the  value  of  r0.  Microsecond  accuracy  in  t„  is  adaquate  for  this  purpose. 
With  t„  known,  the  local  epoch  can  also  be  converted  to  GPS  time.  If  t,  is  the  local  time,  the  corre¬ 
sponding  epoch  converted  to  GPS  would  be  t* 

where 

t*  =  t,  +  t.  -  T„ 

In  this  equation,  the  senses  in  Figure  2  are  used  to  determine  the  signs. 


DOPPLER  OBSERVATIONS 


The  Doppler  observation  consists  of  two  measurements:  the  Doppler  count  is  an  integer  that 
represents  the  accumulated  number  of  Doppler  cycles  and  is  recorded  at  a  regular  interval  7~;  and  the 
time  interval  measurement  is  equivalent  to  a  fractional  cycle  count.  The  interval  measured  is  be¬ 
tween  each  T epoch  and  the  next  positive  going  zero  crossing  of  the  analog  Doppler  signal.  Both  of 
these  measurements  are  illustrated  in  Figure  3. 


The  signal  that  is  counted  is  the  difference  between  an  offset  frequency  /,,  and  the  true  Doppler 
/,,.  Since/,,  -  28750  Hz,  a  count  of  zero  Doppler  would  register  1,725,000  counts.  The  fractional 
count  is  combined  with  the  integer  count  in  the  following  manner.  Referring  to  Figure  3,  the  time  in¬ 
terval  that  represents  the  integer  count  is  T  +  t,  ,  -  t..  Therefore,  the  average  frequency  is  in  ( T 
+  t,-,  -  t , )  Hz.  From  the  average  frequency,  the  approximate  number  of  counts  in  the  interval  T 
can  be  determined.  This  count,  in, ,  is 

s 


mT  =  m  +  (t,-,  -  T,)  _ _ _ 

T  +  T,-,  -  T, 

The  error  introduced  into  mT  by  using  the  average  frequency  is  discussed  in  Appendix  C. 

The  Doppler  data  consist  of  two  simultaneous  observations:  one  from  Li  and  the  other  from 
L2.  These  two  frequencies  will  be  refracted  differently  by  the  ionosphere  because  it  is  a  dispersive 
medium.  Solution  of  the  resulting  simultaneous  equations  allows  a  correction  term  to  be  evaluated, 
and  substitution  of  this  term  into  either  of  the  original  two  equations  will  produce  a  corrected  obser¬ 
vation.  The  mechanics  of  this  process  will  be  developed  in  the  sections  that  follow. 


CLOCK  MODEL 


The  clock  model  used  in  this  development  will  be  expressed  as  a  phase  function  <!>(/). 

<D(r)  =  <)>  +  2n(v0  +  v(f))  t  +  r)(0  (1) 

where 

<t>(r)  =  total  phase  accumulated  since  t  =  0 
<t>  =  phase  at  /  =  0;  i.e.,  0(0)  =  <|> 

Vo  •  =  nominal  oscillator  frequency  (constant) 
v(r)  =  deterministic  frequency  variation 
r)(t)  ~  random  phase  noise  term 

This  form  of  the  phase  will  represent  both  the  local  standard  reference  as  well  as  the  satellite 
reference. 

The  instantaneous  frequency  is 
where 

<t>U)  =  2n(v0  +  v(t ))  +  2  nv(/)t  +  n(/)  (2) 

and  at  r  =  0, 

<t>(0)  =  2n(v„  +  v(0))  +  >)(0) 

The  argument  I  will  represent  GPS  time  as  kept  by  the  MCS.  All  other  clocks  will  have  offset  and  rate  dif¬ 
ferences  when  compared  with  the  reference  standard.  A  clock  keeping  perfect  GPS  time  can  be  represented  by 
a  phase  function  of  the  following  form: 

2  nv„t  (3) 
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where  v0  is  a  known  invariant  constant.  Time  is  kept  by  keeping  track  of  the  phase  <t>GPS(r).  A  time 
interval  would  be  represented  by  the  difference  of  two  phases  at  different  instants. 

h  -  1 1  =  _ J__  -  4>(il>s(M]  seconds 

2nv,, 

With  the  clock  model  defined  as  above,  let  the  satellite  clock  be  represented  by 


<*>,(/)  =  f  +  2n(v„,  +  v.(r ))  /  +  i ^,U)  (4) 

and  the  local  clock  by 

<M t)  =  +  2n(v„, +  v,(t))  t  +  ni(0 


DOPPLER  SIGNAL  MODEL 


The  signals  recovered  will  be  the  transmitted  signals  perturbed  by  the  ionosphere  (troposphere  and 
relativity  effects  will  be  ignored  in  this  development).  The  ionosphere  will  be  introduced  as  a  time  delay 
similar  to  the  geometric  propagation  time.  Let  the  total  propagation  delay  be  represented  by  Tp 

where 

tp  =  _!_  (r  +  I)  and  r(t)  =  slant  range  from  the  satellite  to  receiver 
c 

The  second  term,  /((),  is  the  ionospheric  delay. 


Signals  transmitted  at  t,  are  received  at  a  later  time  tT  +  xp,  or  conversely  the  signals  currently 
received  were  transmitted  earlier  at  tK  -  rp.  Before  transmission.  Equation  4  was  multiplied  by  a 


constant  q,  for  broadcast  on  L, .  Similarly,  Equation  4  was  multiplied  by  for  broadcast  on  L2. 
The  received  signals  are  represented  by  Equations  6a  and  6b. 

<t>„(/„-Tp)  =  <7,[f  +  2n(v„, +  v,)(/„-Tp)  +  r],(/w-Tp)]  (6a) 

Q'lU*  ~  Tp)  =  (h  [f  +  2n(v„,  +  v,)Ur  ~  Tp)  +  r),Un  ~  Tp)]  (6b) 

In  the  receiver,  two  similar  signals  are  derived  from  the  local  clocks;  these  are  represented  by 
Equations  7a  and  7b. 

<*>i  >(/« )  =  ei  I^k  +  2n(v„, +v,W« +  »],(/«)]  (7a) 

4>/i(fw)  =  Pi[+«  +  2rt(v0;  +  v,)/„  +  rj,(rR)]  Ob) 


The  receiver  mixes  the  locally  generated  signals  with  those  received  and  derives  an  intermediate  fre¬ 
quency  which  is  then  processed.  This  intermediate  frequency  phase  function  will  be  written  as  the 
difference  between  Equations  6  and  7: 


<*>A.(/«)  =<M/k  -  Tp)  -  <t>L,VK) 

=  t7,<t>,  -  e,<)>K  +  2n[<7,(v„,  +  v.)(f*  rp)  -Q,(v0,  +  v,  )t„\ 

+  (}xnXtK-rp)-elr)L(tH)  (8a) 

♦«(',.)  =  <Mf*  ~  t„)  -  4>tl(/„) 

=  <72<|),  -  e2<j>*  +  2n[<72(v0,  +  v.)(/«  -  Tp)  -e2(v0;  +  v,  )tH] 

+  qif]s(tK-yll)-Qlr]l(lH)  (8b) 

The  Doppler  observation  is  essentially  a  measurement  of  phase  at  two  different  times.  For  ex¬ 
ample,  from  Equation  8a,  the  phase  difference  T  seconds  apart  would  be 

+  71 -<t>fl(r„) 

A  detailed  expansion  is  written  below  as  A d>: 

A<D  =  <Dm(/«  +  T)-QnitH) 

AO  =  2n\qy\v„,  +  v,(/„  +  71][/„  +  T  -  t p(/„  +  71]  -  qt  [i/„,  +  v, (/„)](/„-  Tp(/„)] 

~  e.[vo/  T  Vi\tH  +  71][r„  +  7']  +  (?2(v<»  T  Vj.((u)]  t n } 

+  <7.  [V'Vk  +  T  ~  Tp(/«  +  71)-r],(7«-Tp(?„))]-  +  7')-rj,(f„)]  (9) 


This  expression  does  not  contain  <7,$.  or  e,$«  because  these  are  constants  and  disappear  in  the  sub¬ 
traction.  In  like  manner,  if  there  were  a  channel  bias  term  in  Equation  8a,  it  would  also  not  ap¬ 
pear  in  Equation  9.  Therefore,  this  observation  of  AO  is  not  sensitive  to  a  channel  bias  as  long  as  the 
bias  remains  constant  over  the  interval  T. 


It  will  be  assumed  that  the  frequency  drifts  contained  in  v,  and  v,  are  small  enough  so  that  the 
differences  over  a  time  T can  be  neglected.  To  demonstrate  this,  suppose  v.(t)  is  on  the  order  of  4  * 
10'"  Hz/s.  This  is  a  representative  number  obtained  from  the  specification  of  the  HP5065A  Rubidium  fre¬ 
quency  standard’.  The  frequency  drift  for  this  unit  is  given  as  less  than  ±  1  *  10  "  parts  per  month.  Multiply¬ 
ing  this  by  10.23  MHz  and  dividing  by  the  number  of  seconds  per  month  produces  the  number  stated  above. 
Over  the  interval  T.  the  difference  <7,  [v,(/n  +  T)(/„  +  71  -  v.(tn)  f*l  equals  4  x  10'"<7,  T.  Using  T  =  60sand 
q,  =  154,  the  resulting  phase  drift  is  2n(4  *  I0'")(I54)(60)  =  2.3  *  lO"  radians,  which  is  small  enough  to  ig¬ 
nore.  However,  the  random  variations  represented  by  17(f)  must  be  considered.  Frequency  instability  is  usual¬ 
ly  represented  as  a  fractional  frequency  variation  A/// over  an  interval  T.  This  implies  a  phase  instability, 
which  can  be  illustrated  if  the  following  operations  are  performed. 


Let 

V  <71  =  A 

/ 

Then 
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4/<  T)  =  fk 


and 

/ 

dt  n  ¥  ( T ) 

/ 

Therefore,  ihe  drift  due  to  the  random  contributions  from  the  local  reference  frequency  standard 
(assuming  T  ■=  60s,  A  /./<D  -  I  *  10  './  -  10‘H/  and  e,  -  315)  is  2n(10  u)  (315)  <5  *  H)*) 

(60)  =  0.47  rad.  This  magnitude  of  phase  instability  is  considerable  and  constitutes  a  major  source 
of  error  in  determining  the  phase  difference  A4>.  This  phase  noise  will  be  collected  and  represented 
by  the  parameter  ip  where 

4>  =  r/,h.u«  t-  r-T  {,(„  +  7~>)  -  rj.  (/„  —  T^(/« )))  -  t>!  (Hj  (/«  +  H  -  q, 0* )|  (10) 

Since  v.(r)  and  v,U)  were  shown  to  vary  little  in  an  interval  T.  equation  9  can  be  simplified  to 

A<J>(/k  +  T)  -  2n|</l(v„.  +  v.) \T  *  r  (/„)  jjl„  ♦  D)  -  e,(i/,„  +  v, )  T\  +  tp 

Division  of  A <t>  by  2n  produces  a  number  V.  which  is  the  Doppler  count  (th-‘  number  of  cycles  in  the 
interval  T). 


T)  -  r\(iH )  =  f  A/r/r  fk  J 


,V,  =  A4»0k  r  T)  r/,(t/„.  +  v.)|7'  ♦  t,(u)  tJim  ♦  7'))  -  c, (w,„  +  v, )  T  +  ^  (11) 

2n 

This  is  easily  rearranged  to  place  the  difference  in  r  \  on  the  left  hand  side. 


zn 


t^Ok)  -  T,(/„  ♦  7)  1  f.\,  -  (</,(v„.  +  v,)T- e,(v«,  +  V,  )D  -  V 

</■(*'.*  +  v.)L  2n 


(12) 


In  Equation  12,  V,  is  the  observation.  The  rest  of  the  terms  in  the  brackets  represent  that  portion  of 
the  count  that  is  due  to  the  Doppler  offset  and  the  instability  of  the  frequency  reference. 


I  el  us  hypothesi/e  that  the  propagation  delay  r .(/)  is  composed  of  the  geometric  slant  range 
r(r)  and  the  ionospheric  delay  Hi),  which  were  introduced  earlier.  If  this  is  substituted  into  Equation 
12,  the  result  is 


i  {I  h  )  ~  r*(  /  h  *  7)  *  / 1  ( r « )  /  Hi  h  *  7)  c 
In  like  manner,  an  expression  for  /  ,  can  be  developed: 
r(t 


V,  (</,(!,„  +  V  )  7  «■  V,  |T|  \p,  2 IT 

<{Av„.  +  v  ) 


L  (1 3a) 


t„)  rV„  *  7)  *  lAt„)  /,(/„  *  7)  «J-V  +  v.)  f  qAv.„  +  „,)7]  -  \p,  In 

\  I/;(V„.  +  V,  ) 


(13b) 


Since  the  geometric  range  r(/)  will  be  common  to  the  two  frequencies,  it  can  be  eliminated  bv  sub¬ 
tracting  Equation  13a  from  Equation  13b. 


/,(r«)  /,Uh  *  7  )  )/,</„)  -  /,(/«  ♦  7  )j 


'.V, 

V, 

*  (v,„  +  v,  )7 

t?. 

-  Qi 

1 

(h>< 

(14) 

.  <(< 

<h 

<7. 

<h 

2  n 

V/. 

•t>j\ 
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This  expression  equates  the  part  of  the  range  that  is  frequency  dependent,  on  the  left,  to  the  obser¬ 
vations  N, .  N,,  and  T on  the  right.  The  ionosphere  can  be  modeled  by  expressing  the  quantity  /(/) 
by  an  integral  expression'. 


/</«)  = 


-1  f  u /(r)  rdr 

2wl  Jv»>  V^:k1)'n 


(15) 


where 


/•,(/„)  =  satellite  position  at  /« 
rjr„)  =  recc. . cr  position  at  /*  +  r. 
u),  =  2iu/,(v„.  +  v.) 


A  similar  expression  can  be  written  for  /,(/«  +  71  having  the  limits  r,(/K  +  71  and  rt(tK  +  T).  The 
transmitted  frequency  appears  only  in  the  constant  coefficient  u>, .  The  difference  /,(/*)  -  /,(/«  +  T) 
can  be  written  as  the  difference  of  two  integrals  like  Equation  (15). 


M/»)  -  I<Uh  +  71  - 


oj  ‘(r)  rdr 
(r^JJTJTTT 


oj,‘(r)  rdr 
(rl  -  k1)'1 2 


(16) 


The  expression  in  braces  will  be  identical  for  l and  so  it  can  be  redefined  as  a  constant  /  such  that 


.  r  c u‘{r)rdr  f  wf(r)  rdr 

(r‘  A1)''7  J 

Then  Equation  16  can  be  simplified  to  the  following  form: 


/.(/-« )  l,{rH  *■  71 
1  ikewise, 

lAr„)  I t(r h  »  71 


(17) 


Inserting  these  into  Equation  14  and  striving  for  /  produces  an  empirical  result  for  the  value  of  the 
two  integrals  defined  in  Equation  1 7, 


(v„.  -  V  )  J 

f  A  1  -  A;  +  +  l/(  )  J 

Vi  -  Ql  ’ 

</;  <tl 

I 

\  '/.  <h 

//.  <h  . 

1 

Vi  4^1 

2rr 

//■  (h\ 

(18) 


This  can  now  be  substituted  into  either  Equation  13a  or  Equation  13b  to  produce  an  ionosphere- 
corrected  observation  of  the  range  difference  r(/„)  -  r(iH  +  71.  These  are  written  as  Equations  19a 
and  19b: 


V,  -  |r/,(i/„,  +  v  )J  -  +  v,  )  1 ) 

i/Av„.  +  v  ) 


ip,  2n 


+ 


r[tH)  r{iH  *  7) 


/ 

Kn\/,'(v„.  +  v  )' 
10 


(19a) 


(19a) 


r(tK)-r(tH  +  7)  =  r  J  Ni~  l<h(v «'  +  »-)r-  e2(vot  +  vt)71  -  y2/2 n 
1  <72(v0>  +  v.) 


8n2gi(v„,  +  v,)2 


ERROR  ANALYSIS 


Now  that  the  conversion  from  raw  Doppler  counts  to  range  difference  has  been  developed,  it 
will  be  informative  to  investigate  the  sensitivity  of  these  expressions  to  errors  in  the  measurements. 
The  last  term  in  Equation  19  can  be  combined  with  Equation  18  to  eliminate  some  of  the  constants. 
This  new  parameter  is  defined  in  Equations  20a  and  20b: 


J„ 


Jh 


I  =  eg  j  W.  _  Nj 

8tt2<7,2(v„,  +  v,)2  (q{  -  r/22)(i/„.  +  I/,)  q2 


T 

[£a 

-  £i 

1  -  — - 

[>> 

_ 

Qy 

J  ^ 

U 

9*. 

I 

= 

aft 

Nj 

8n2qi(v0.  +  v,)2 

u/t  - 

qi  )(v„,  +  v.) 

l  Q' 

Qi 

T 

-  £il 

1 

‘Vi 

_ 

<h 

2n 

<7. 

J. 

+  (v0i  +  v, ) 


+  (l/„l  +  v, ) 


(20a) 


(20b) 


The  parameters  that  will  be  considered  to  be  uncertain  are  v..  v, ,  /V, ,  Nt.  q>,.  qi2  and  T. 


Special  consideration  must  be  given  to  qj,  and  qj2 ,  since  they  are  highly  correlated.  The  form  of 
Equation  10  can  be  used  to  describe  both  if,  and  q>2. 

4*.  =  r/,(r).(r„  +  T- rf>)  -  -  r,)|  -  c , |rj, (/«  +  7~)  -  rj, )| 

=  '//[*?.(/«  +  T -  Tt  )  p,(/w  -  7(,)|  -  €>/!>?/(/*  +  T)  -  r)A> »)] 

Using  these,  the  quantity  ,  </,  -  \p2  </2,  which  is  common  to  both  J„  and  J,,,  can  be  written: 

9'.'  -  lf2i  =  \rjAln  +  T- t/;)  -  r\,{lH  -  t(i)|  -  e'  1*1  <(t«  +  T)  -  *7, (/„)] 

Q\  <h  (h 

-({*?,</«  +  T -  t,)  -  r\\tK  -  T„))  -  Ci  |*),(f«  +  T)  -  *i,(/k)H 

<h 

Here  the  satellite  phase  noise  contribution  combines  to  add  to  zero,  which  leaves  only  a  reduced 
contribution  from  the  local  frequency  reference. 


II 


Vl  _  V2  = 
Q>  Qi 


(e±  -  M 

\Qi  Q>) 


(ni(tn  +  T)  -  ndtR)) 


=  -1.06  x  10'5[rjt(/«  +  T)  -  rjiUn)] 

Thus,  in  addition  to  partials  of  Ju  and  Jh  with  respect  to  u>t  and  ip2,  the  partials  dxpjdr,  and 
dy2/dr,  will  be  defined. 


dip, 

=  q,  /£i 

dr, 

U* 

dy2 

=  -  3t v, 

<h 

<7. 


3  r,  dr, 

where 

1  =  hi(tR  +  T)  -  rj,(rR)] 

The  partials  of  J  with  respect  to  the  seven  parameters  are  written  below: 


dj„ 

— 

-cql 

l 

-  ^ 

dv% 

li  - 

(V 

>.  + 

V.)2 

_ 

{<?• 

92 

T 

ii  - 

o£ 

1 

_  vT 

q> 

q\ 

2n 

U- 

dJh 

_ 

-cqi 

I 

f". 

Ni 

dv, 

i  -  qi 

(v0,  + 

v.)2 

1?. 

q2 

T 

[£i  - 

Qj_ 

1 

- 

k* 

qi 

2n  [7, 

qi 

d  Ja 

cql 

T 

e.  _  £2' 

3v, 

(qi  - 

ql)(v0. 

Tv 

) 

<7l  <?2 

dJa 

= 

cqi 

T 

-  £2' 

dv, 

(qi  - 

ql)(vn< 

+  V, 

) 

//.  q  I* 

dja 

_ 

cql 

9N, 

(qi  - 

ql)(vlu 

+  V 

)q> 

d  jh 

_ 

cq, 

a/v, 

(qi  - 

ql  )(v„. 

~TT 

) 

dJ„ 

_ 

-cq>_ 

d  n2 

(qi  - 

Tv, 

)~ 

i  +  (vn,  +  v, ) 
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dJh 

= 

3Ni 

(9?  -  <7i)(v<>,  +  v,)<72 

37u 

=  -cql 

3i^i 

(q<  -  ql)(v o,  +  v,)2nqt 

3  Jb 

=  -cq, 

3tp, 

(<7?  -  <7il)(v0,  +  v,)2tt 

37„ 

=  cq2 

3t^i 

(<7?  -  qi)(vo,  +  v,)2n 

3  A 

=  c<72 

34»j  iq>  -  Ql)(vo,  +  v,)2nq2 


dJ°  =  cqj{.-0l  +  vL ) _  [c,  _  p/j 

dT  (qi  -  7?Kvo.  +  v.)  [<7,  q2\ 

dJ*  = _ <y,;(vot  +  Vl>  [e.  -  £j"| 

dT  (<??  -  v,  ti.  v,  +  v,)  [ q ,  q2\ 

The  error  in  J  due  to  these  parameters  is  obtained  from  the  total  derivative: 


dj  =  AL  dv ,  +  AL  dvt  +  AL  dN  +  AL  rf^+  AL  dT  (2i> 

3v,  3v,  dN  dip  dT 

If  values  are  assigned  to  the  constants,  the  sensitivity  of  J  to  each  of  the  parameters  can  be 
evaluated.  The  appropriate  constants  are 

Vo.  =  10.23  x  10*  Hz 
v,  =  0 

Voi  =  5.00  x  10*  Hz 
v,  =  0 

c  =  2.99792458  x  10'  m/s 
q,  =  154,  </,  =  120 
q,  =  315.07825,  Gj  =  245.51425 
T=  60  s 

N,  -  Nj  =  1725000 
=  0 

Using  these,  the  numerical  values  for  the  partials  are 

AL  =  -  0.0992  m/Hz,  AA  =  -  0.1634  m/Hz 
3v,  3v, 

AA  =  0.0288  m/Hz,  AA  =  0.0474  m/Hz 
3v,  3v, 
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=  0.2941  m/count,  d-A.  =  0.4844  m/count 
ON,  ON, 


=  -  0.3775  m/count,  £ =  -  0.6217  m/count 

on2  on, 

=  -  0.0468  m/rad,  d-A  =  -  0.0771  m/rad 
0  i  Oip, 

a-/-  =  0.0601  m/rad,  _££i  =  0.0989  m/rad 

0yj2  0t//2 

.££l  =  2396.0128  m/s,  _£A  =  3946.1000  m/s 

or  or 

In  a  similar  fashion,  Equation  19  allows  the  definition  of  A  and  A  as  the  uncorrected  A,  and 
A2  range  difference  observations: 


e 

II 

-A. _ 

In,  -  [<7,(v0. 

+  v,)r-  e,(v0/.  +  v,. 

)  r]  - 

Vi/2n] 

l 

(22a) 

1 

l 

<7.(v«,  +  v,) 

J 

A  =  c  J 

|n2  -  [<72(v0. 

+  v,)T  -  e2(v0/.  +  v. 

)71- 

t/>2/2n] 

• 

(22b) 

1 

l 

q  2  ( v„.,  +  v,) 

J 

partials  of  D„  and  D , 

>,  with  respect  to  the 

same 

set  of  parameters 

are 

0A  s 

=  -cr 

_  c 

In, 

-  1  a,(vn.  +  vAT 

-  e.(v0/. 

+  v2)r]  - 

V.  1 

3v, 

Vo.,  +  V., 

<7 1  ( v0,  +  v,)2 

1 

2n  J 

a 

=  -cT 

_  c 

Jn2 

-  fO,(Vn.  +  vAT 

-  e2(v0I 

+  vjr]  - 

Vi  1 

3v, 

V«,  +  V, 

q  2(v0,  +  v,)2 

2n  j 

0  A 

II 

i 

<8 

“i 

Ov, 

9'(Vo,  +  v,) 

0  A 

=  -cqz  r 

Ov, 

q  2  ( v0,  +  v,) 

0  A 

=  0A  =0 

ON, 

0  n2 

OA 

=  c 

OA  =  c 

ON, 

0N2  <72(v„,  +  v.) 

0D„ 

—  -c 

Oif», 

2nq,(v„,  +  v,) 
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dDb  = 

—c 

dtp2 

2n<72(v0,  +  v5) 

CD 

& 

» 

li 

dDa  =  o 

dtp, 

dt pi 

dD„  = 

c  [-  <7i(vo.  +  v.)  +  G,(v0i 

+  Vi)] 

dT 

Q  i(v0,  +  v.) 

dDh  = 

C  I"  <72(Vo,  +  v.)  +  Qi(.Vol 

+  Vi)) 

dT 

Qi(v0,  +  v() 

The  sensitivity  of  D  to  each  of  the  parameters  can  be  evaluated  using  the  same  set  of  constants 
as  before.  The  partials  with  respect  to  v ,  are  a  function  of  the  actual  Doppler  count  and  are  a 
minimum  when  the  Doppler  frequency  is  zero.  This  is  because  the  terms 

<7i(v0,  +  v,)  T  -  e,(v0i  +  vl)T 

and 

<7i(v0,  +  v,)  T  -  Gi(v0l  +  v,.)T 

equal  1,725,000.  When  N,  =  Ni  =  1,725,000,  the  term  in  braces  is  zero.  In  order  to  represent  the 
worst  case,  /V  will  be  adjusted  to  the  equivalent  of  a  ±4000  Hz  Doppler  frequency  in  these  two  par¬ 
tial  only. 


dDa 

dv, 

=  -1758.3036  m/Hz, 

dDh  = -1758.3048  m/Hz 
dv. 

d  Da 
dvL 

=  -3597.4143  m/Hz, 

dD*  =  -3597.3957  m/Hz 
dvL 

dDa 
d  N, 

=  0.1903  m/count, 

d  Dh  =  0.2442  m/count 
dN2 

dDa 
dtp , 

=  0.0303  m/rad. 

dD>  =  0.0389  m/rad 

dtp  2 

dDa 

=  -5470.9  m/s. 

dD "  =-7021.0  m/s 

dT 


dT 


A  comparison  of  the  values  of  the  partials  for  D  and  J  indicate  that,  in  all  cases,  the  greatest  con¬ 
tribution  to  the  observation  of  a  range  difference  will  come  from  the  D's  not  the  f  s. 


Equations  19a  and  19b  are  equivalent;  therefore,  the  range  difference  observation  can  be  ob¬ 
tained  from  either.  However,  inspection  of  Equations  20a  and  20b  shows  that,  for  a  given  error  in 
any  of  the  quantities  in  the  braces,  the  resulting  error  in  J„  will  be  less  than  in  Jh.  This  is  because  of 
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the  different  multiplying  factors  involved.  Also,  Equations  22a  and  22b  indicate  that  an  error  in  the 
count  N,  time  interval  T,  or  in  the  phase  noise  ip.  will  be  less  in  Da  than  in  Db.  Again,  this  is  because 
of  the  multiplying  factors.  This  conclusion  is  also  apparent  in  the  values  of  the  partials  listed  above. 
In  all  cases,  the  partials  of  Da  and  Ja  are  less,  respectively  than  those  of  Db  and  Jb.  However,  the  er¬ 
ror  propagated  into  the  ionospheric  corrected  range  difference  observation  due  to  errors  in  Not  ip  is 
the  same  irrespective  of  whether  Da  or  Db  is  used  as  the  basis  of  the  measurement. 

This  can  be  illustrated  in  the  case  of  an  error  in  N  by  employing  the  appropriate  partials.  The 
partials  derived  from  Equations  19a  and  20a  are 


Ar  = 

dD“  AN,  + 

dJa  AN,  + 

&  AN* 

dN, 

dN, 

dN2 

Ar  = 

c 

AN,  + 

cql 

q,(v0,  +  v.) 

(qf 

-  ql)(v0,  +  vs)q, 

Ar  =  _ £ _ AA\  +  _ _ AN,  +  _ £  cQi) _ AN 2 

q i(v0. +  v«)  (q1,  -  ql)(vB>  +  vs)q,  (q\  -  q\)(v0s  +  v,) 

Assuming  that  AN,  =  A N2,  this  results  in 

Ar  =  c  f  1  AN 

[(<7?  -  <7f  )(v0v  +  v.)  J 

A  similar  expression  can  be  written  from  Equations  19b  and  20b: 

Ar  =  JiJ  AN2  +  lit  AN,  +  JA  A N> 
dN2  dN,  d  n2 

Ar  =  AN2  +  c<?» _ AN,  -  _ ££i _ AN2 

q2(v„ ,  +  v,)  (qi  -  ql)(v0,  +  v,)  (q1,  -  q})(v 0s  +  v.) 

Assuming  that  AN,  =  A N2,  this  also  results  in 

Ar  -  c  Si  ~  Qi  AN 

.  (q f  -  <72)(v„,  +  v,)  _ 

A  similar  result  is  obtained  when  the  error  parameter  is  ip.  Consequently,  if  the  primary  error  source  is  due  , 
to  uncertainties  in  N  or  ip,  either  expression  produces  an  equivalent  error  in  the  final  range  difference 
calculation.  Since  this  is  the  case,  the  next  section  will  dwell  on  evaluating  the  errors  from  D„  and  Ja. 


EXPECTED  OBSERVATIONAL  ERROR 


An  estimate  of  the  error  in  each  observation  can  be  made  by  using  the  preceding  partial 
derivatives  and  the  expected  errors  in  each  of  the  parameters  considered.  These  expected  errors  will 
be  based  upon  the  considerations  discussed  in  the  following  paragraphs. 


The  frequency  offset  v,  will  be  expected  to  drift  at  a  rate  of  about  4  x  ]0’18  parts/s.  Using  the 
nominal  frequency  for  v„,,  this  is  about  4.1  x  10~"  Hz/s,  or  2.5  x  IO’9  Hz  in  a  60-s  observation 
period.  The  frequency  offset  v,  will  be  expected  to  vary  less  than  1.2  x  10~18  parts  of  v0/.,  or  3.5  x 
IO"10  Hz.  This  latter  number  is  from  the  HP5061A-004  specification. 


The  Doppler  count  accuracy  is  limited  by  the  resolution  of  the  fractional  time  interval  counter. 
At  the  nominal  offset  frequency  of  28750  Hz,  the  ±  10-ns  resolution  of  the  time  interval  counter  is 
equivalent  to  approximately  2  x  10’8  x  28750  of  a  count.  This  is  5.75  x  10~4  cycle,  or  3.6  x  IO’3  rad. 
The  white  noise  from  the  receiver  is  on  the  order  of  0.1  rad  rms,  so  the  limiting  factor  is  not  the 
counter,  but  this  receiver  noise.  Thus,  the  Doppler  count  error  will  be  taken  as  0.1  rad,  or  0.016 
count,  in  the  rest  of  this  development.  This  error  source  should  decrease  with  increasing  signal 
strength. 

The  phase  noise  ip  comes  primarily  from  the  frequency  standards.  Equation  10  shows  that  it  is  a 
composite  of  contributions  from  the  satellite  and  local  references.  The  satellite  will  be  assigned  an 
Allen  variance  over  a  minute  of  2  x  IO’12  (Reference  4),  while  the  local  reference  will  be  assigned  1  x 
IO’12.  The  phase  variation  over  1  min  can  be  estimated  from  Equation  10  in  the  same  manner  as 
before.  This  expression  when  evaluated  becomes 

+  T)  -  »»,(/«)]  =  2n(2  x  IO'12)  (154)  (10.23  x  106)  (60)  =  1.19  rad 
In  a  similar  fashion,  the  result  for  the  local  clock  is 

e,M'«  +  T)  -  n,(t«)]  =  2n(l  x  10"2)  (3.15.08)  (5  x  106)  (60)  =  0.59  rad 
From  both  of  these,  the  magnitude  of  Ai+i  =  1.33  rad. 


j 

1  4 


The  accuracy  of  the  Doppler  count  interval  is  correlated  with  the  count  itself,  since  the  interval 
is  triggered  to  begin  and  end  by  the  zero  crossings  of  the  Doppler  count.  Since  this  signal  has  noise 
attached  to  it,  there  will  be  a  corresponding  uncertainty  in  the  duration  of  the  interval.  Using  the 
0.1-rad  figure  as  representative  of  the  noise  and  the  duration  of  the  nominal  Doppler  cycle  as  3.5  x 
IO’5  s,  the  time  corresponding  to  0.1  rad  of  this  frequency  is  5.5  x  10’7  s. 

Finally  the  error  in  the  observation  will  be  calculated  using  the  total  derivative  as  in  Equation 

21. 


AJ„  =  Av,+  Av,  +  JLL  AN,  +  A N2  +  ALl  A ip,+  Hjl  Aip2  + 

3v,  3v,  3N,  3  N2  dip,  d\p2 


3  J. 


AT 


3  T 


The  magnitude  of  A J„  is  obtained  from  the  square  root  of  the  sum  of  the  squares  of  each  of  the 
terms  above.  A  similar  expression  will  be  used  for  A D„. 


A  D„  =  dD“  Av,  +  dD-  Av,  +  dD-  AN,  +  dD »  Aip,  +  d  D°  AT 
3v,  3v,  3 N,  dip,  3  T 


The  results  of  the  calculation,  term  by  term,  follow: 


J 

dJ„ 

Av, 

=  -2.5  x  10""  m. 

Av,  = 

-4.4  x  IO’8  m 

3  v, 

d  v. 

% 

3./,, 

Av, 

=  1.0  x  10  "  m. 

dD„ 

Av,  = 

-1.3  x  IO’8  m 

3  v, 

3  v, 

■"  /'  * :  t.  ^ 

imm:  a.-*  ^ 
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a./., 

A.Y, 

0.004'  m. 

^  AN 

0.0030 

111 

a  a, 

a  \ , 

dJL, 

A\, 

-  -0.0060  m. 

•”>  A  \ 

0 

d  Nt 

a  v. 

a./., 

A  40 

=  a  ./„  a  4 1 1 

At]  =  -8  x  |() 

'  Ill, 

d  l)‘  Aip,  =  0.040  111 

dig, 

<54 0  dr] 

dip, 

a  j„ 

A  40 

=  a./„  dip. 

Aq  =  -8  x  10 

'  Ill, 

11 

£ 

< 

** 

dip. 

dip,  dr/ 

dip , 

dJ„ 

AT 

=  0.0013  m. 

AT  = 

-0.0030  1 

m 

3T 

dT 

It  is  interesting  to  see  that  the  aecuracy  of  the  ionospheric  correction  is  limited  not  by  phase  noise 
from  the  frequency  standards,  but  by  the  errors  in  determining  the  Doppler  counts  N,  and  N, .  This 
noise  comes  primarily  from  random  phase  fluctuations  in  the  receiver.  If  N,  and  /V2  are  considered  to  be 
independent,  then  the  magnitude  of  the  noise  expected  on  the  ionospheric  correction  observation  is 


AJ„  = 

f  (U  A.Y,  \ 

■’ +  (3J,  a,Y2  ) 

2  ■ 

\a.v,  / 

\  a  ,y2  / 

. 

AJ„  =  0.0076  m 

The  contributions  from  the  remaining  parameters  are  small  enough  to  ignore. 


The  error  in  the  range  difference  observation,  4.02  cm,  is  dominated  by  the  frequency  standard 
phase  noise.  If  this  error  source  were  reduced  by  a  factor  of  20  or  more,  then  the  limits  would  be  set 
by  the  error  in  determination  of  the  count  and  the  exact  interval  of  the  count.  This  factor  of  20  is 
equivalent  to  improving  the  Allen  variance  from  10  12  to  5  x  10  14  oxer  the  interval  T and  would  be 
comparable  to  the  performance  expected  from  a  Hydrogen  Maser. 


The  levels  of  error  in  D „  and  J„  predicted  by  the  preceding  analysis  have  been  empirically 
verified  by  actual  data.  NAVSTARS  1  through  4  (SV’s  •»,  7,  6,  and  8)  were  operational  while  the 
NGRS  was  at  Yuma  in  1979.  The  Doppler  data  that  were  obtained  from  all  four  satellites  were  pro¬ 
cessed  so  that  20-min  segments  were  fitted  to  a  third-degree  polynomial  using  least  squares.  The 
residuals  of  this  fit  were  used  to  determine  an  rms  that  is  considered  to  be  representative  of  the  ran¬ 
dom  fluctuations  of  the  Doppler  observations.  Uncorrected  L,  Doppler  data  corresponding  to  £>,, 
and  the  ionospheric  correction  to  L  , ,  corresponding  to  J, ,  were  processed  in  the  same  way.  The  rms 
values  of  the  residuals  for  D  are  plotted  vs  date  in  Figures  4  through  7,  and  the  residuals  for  J„  are 
plotted  in  Figures  8  through  II. 


Throughout  these  44  days,  the  local  frequency  reference  behaxed  normally,  and  so  the  jumps  in 
the  calculated  rms  (e.g..  those  in  Figure  4  between  days  45  and  54)  are  due  to  misbehavior  of  the  SV4 
frequency  standard.  Note,  however,  that  during  these  nine  days  there  is  no  evidence  of  anything 
amiss  in  Figure  8.  I  his  illustrates  that  the  ionospheric  correction  is  not  sensitive  to  the  phase  noise  of 
the  Irequencx  standards.  This  was  the  result  predicted  by  the  analysis  of./,,.  An  especially  graphic 
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figure  5.  2()-Min  /  ,  Range  Difference  RMS  (cm)  Vs  Date. 
GPS  SV-6,  Yuma  Proving  Grounds.  Days  36-80  1979 
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Figure  7.  20-Min  Range  Difference  RMS  (cm)  Vs  Date, 
GPS  SV-8,  Yuma  Proving  Grounds,  Days  36-80  1979 
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Figure  8.  20-Min  Ionospheric  Correction  to  L,  RMS  (cm)  Vs  Date, 
GPS  SV-4,  Yuma  Proving  Grounds,  Days  36-80  1979 
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Figure  10.  20-Min  Ionospheric  Correction  to  L,  RMS  (cm)  Vs  Date, 
GPS  SV-7,  Yuma  Proving  Grounds,  Days  36-80  1979 
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Figure  II.  20-Miti  Ionospheric  Correction  lo  /.,  RMS  (cm)  Vs  Date, 
CPS  SV-8.  Yuma  Proving  Grounds,  Days  .16-80  1979 
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demonstration  of  this  is  available  from  SV7  in  f  igures  6  and  10.  The  rms  of  the  residuals  of  D„  for 
this  satellite  were  poor  throughout  the  44-day  span;  however,  the  ionospheric  correction  residuals  in 
Figure  10  are  unaffected  and  are  comparable  with  those  from  any  of  the  other  satellites.  This  feature 
is  certainly  a  plus  for  anyone  interested  in  performing  ionospheric  measurements,  because  it  means 
that  they  need  not  invest  in  a  very  low -noise,  and  consequently  expensive,  local-frequency  standard. 

Since  the  error  in  D.,  is  primarily  due  to  At p  and  the  error  in  J„  is  due  to  A/V.  they  should  be  in¬ 
dependent.  Therefore,  the  predicted  error  in  range  difference  would  be  equal  to 

A/-  =  [(AD,.)'  +  (A7„)') 

Using  the  numbers  derived  previously,  the  range  difference  observation  error  should  be  about  4.0 
cm.  This  is  in  excess  of  what  is  observed  for  satellites  4,  6,  and  8,  but  less  than  that  for  satellite  7.  If 
the  primary  error  source  is  due  to  the  random  phase  fluctuations  of  the  frequency  reference,  then 
assuming  that  the  phase  fluctuations  of  the  local  reference  are  represented  by  o, ,  the  estimates  of 
the  satellite  values  can  be  made  using  the  method  previously  outlined.  The  numerical  value  for  o, 
will  be  taken  from  the  Allen  variance  at  T  ~  60  s  (Table  2).  With  this  value  fixed,  a  number 
representative  of  the  satellite  frequency  reference  can  be  calculated  using  the  following  relations: 

A/;„  =  dl) ■■  AVl 
dyj, 

and 

At^i  =  {r/fh.(/w  *  t  )  m (/«)!“  +  *  T)  )l  ‘ f ' 

for  this  application 

AV,  -  2n  x  60  x  |0-  ((10.2.1  x  154  x  4)  +  ,s  x  115.07825)  | '  '  a, 

Solution  for  £  and  substitution  for  the  average  A/>  obtained  from  Figures  4  through  1 1  over  the  en¬ 
tire  44  days  gives  the  results  listed  in  Fable  2. 


Table  2.  Satellite  Phase  Noise  Derived  From  Doppler  Observations 


sv 

A/>, 

Satellite  Phase  Noise 

Assumed  1  ocal 

(cm) 

£o, 

o, 

4 

2.^8 

1.14  x  10  '' 

1  X  10  " 

6 

2.57 

1 .04  x  10  ' 

1  x  10  ' 

7 

10.61 

5.81  x  |0 

1  X  10  '■ 

8 

1.76 

1.85  *  10  w 

1  X  10  w 

CONCLUSIONS 


The  formulation  for  converting  Doppler  counts  into  range  difference  observations  has  been 
developed  for  the  case  of  the  NGRS.  The  two  frequency  data  were  combined  to  generate  the  first- 
order  ionospheric  term.  When  this  term  is  subtracted  from  the  uncorrected  L,  range  difference,  the 
result  is  the  ionospheric-corrected  range  difference. 

The  sensitivity  of  the  observation  equation  to  errors  in  the  observed  parameters  was  in¬ 
vestigated  by  evaluating  the  partial  derivatives  and  using  them  in  the  expression  for  the  total 
derivative.  Errors  in  the  parameters  were  estimated  and  used  to  determine  the  error  expected  in  the 
range  difference.  These  estimates  were  then  compared  with  examples  of  real  observations  obtained 
when  the  NGRS  was  on  site  at  the  Yuma  Proving  Ground.  Agreement  was  reasonable  and  tends  to 
confirm  that  the  analysis  is  correct.  Further  tests  are  planned  based  upon  the  predictions  that  can  be 
made  from  results  developed  in  this  report. 
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APPENDIX  A 

Frequency  Spectrum  of  the  P  Code  Modulation 
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Define  the  phase  argument  of  a  cosine  function  to  be  9(t)  where 
9(t)  =  oo,t  +  d„  +  A\?(r) 


The  parameters  o>...  0„,  and  k  are  constants,  and  g(t)  is  a  random  gate  function  of  unit  amplitude. 
The  random  character  of  £(/)  is  best  demonstrated  by  Figure  A-l.  A  mathematical  expression  for 
g(t)  can  be  devised  by  the  time  translation  of  an  elementary  gate  function  centered  on  zero  as  in 
Figure  A-2.  Using  this  concept,  we  can  express  the  random  gate  function  as  the  sum  of  translated 
g'(f)' s  multiplied  by  a  random  variable  a,. 

g(t)=  a,g(t-h)  (A- 1 ) 


The  value  of  a ,  is  either  one  or  zero  depending  upon  the  outcome  of  some  random  process  related  to 
the  subscript  /. 
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Figure  A-l.  Random  Gate  Function 
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Figure  A-2.  Elementary  Gate  Function  Centered  On  Zero 


If  we  choose  0„  =  0  and  A 
9(t)  =  (v,t  +  ng(t) 

i 


n.  then  9(1)  becomes 


.11 


(A-2) 


x  utib  bcJUiK-NOi1 


Now  according  to  the  definition  of  g(r),  the  value  of  6(1)  can  be  only  co./  +  n  or  co,/.  This  is  the 
biphase  part  of  the  argument.  When  6(/)  is  introduced  into  a  cosine  function,  the  result  is 


cos  6(1)  =  cos|co.  I  +  ng(/)] 

which  when  expanded  becomes 

cos  6(1)  =  cos  cos  tt g(/)  -  sin  c o.  /  sin  rrgf/) 

Now,  since  ng ( / )  is  either  0  or  n,  the  second  term  vanishes,  which  leaves  just 

cos  6(i)  =  cos  co.  /  cos  rtg(/)  (A-3) 

The  second  cosine  factor  here  is  either  +  1  or  -1  depending  on  whether  the  argument  is  0  or  n.  which 
in  turn  depends  upon  whether  a,  is  1  or  0.  Since  this  is  the  case,  we  can  construct  a  new  factor  to 
replace 

cos  ng(/) 

which  behaves  in  the  same  manner.  Let  this  be  m(t)  where 

/«(/)  =  2 :g(/)  -  1. 

Equation  A-3  written  using  this  new'  factor  is 

cos  6(1)  =  (2 g(/)  -  I)  cos  co.  t.  (A-4) 

The  random  function  m(i)  represents  the  pseudo-random  code  modulating  cos  co.  I.  The  Fourier 
transform  of  the  elementary  gate  function  is 

F{g'(/)1  =  t  _|jn  ^1/2  =  G(co)  (A-5) 

cot/2 

Employing  the  time  shift  property, 

/•)/;(/  -  D 1  =  //(co)  c>  ' 

allows  f  '{a  (I  -  m ) \  to  be  evaluated  from  /-{g '(/)}  as  shown  below  : 


/  jg  (/  -  //t  ) }  =  t  MM  aiT/  “  c’  G{ co)  e 

cot/2 

It  follows  then  that 

/  {g(0  (  =  T  S'n  a,T  ^  X  C/,(*  '  =  (»(co)  X  (/,<' 

cot/2 


( A-6) 


with  t he  implication  i Ik  t  the  contribution  from  any  term  of  the  sum  with  a,  =  0  is  also  zero.  Thus, 
the  frequency  spectrum  of .«(/)  is  a  summation  of  the  form  shown  above  with  an  arbitrary  number 
of  possible  contributors  missing.  The  terms  missing  are  specified  by  the  outcome  of  the  random  pro¬ 
cess  controlling  a  . 


The  basic  frequency  spectrum  (Equation  A -5)  behaves  as  a  sin  x/x  function  multiplied  by  a 
series  of  complex  function  given  by  the  summation  in  Equation  A-6.  The  first  nulls  of G(u')  arc 
when 


n . 


This  frequency  is  10.23  MHz  for  the  GPS  Pcode.  Thus,  the  bandwidth  between  first  nulls  is  20.46 
MHz  and  contains  the  majority  of  the  transmitted  energy.  Since  energy  per  unit  bandwidth  is  pro¬ 
portional  to  | G(o>) |\  a  comparison  can  be  made  by  integrating  |G(o;)|2  between  zero  and  infinity 
and  again  between  zero  and  the  first  null.* 

r  sin  a  dx  =  n  =  j  >->071) 

—x,  "  s 

whereas 

f  sin'vf/iv  =1.41815 

o  *  > 

V 

This  ratio  shows  that  the  fraction  of  energy  between  first  nulls  is  90  percent  of  the  total. 


It  is  interesting  to  see  how  the  gate  function  pictured  in  Figure  A-2  changes  when  G(ce)  is  cut  off 
due  to  the  finite  bandwidth  of  the  transmission  and  receiving  systems.  We  can  approximate  this  ef¬ 
fect  by  defining  a  gate  function  in  the  frequency  domain  fT(u>)  such  that  it  is  unity  between  ±  a>,  and 
zero  elsewhere.  This  function  (pictured  in  Figure  A-3)  is  then  multiplied  by  G(o>),  which  gives  the 
product  G(ce)  H  (ce).  Since  multiplication  in  one  domain  is  equivalent  to  convolution  in  the  other  do¬ 
main,  the  resulting  function  in  the  time  domain  is  the  pulse 

p(G  =  J  n-(r  l  >  e  (V  dl 
Since 


»•(/)=  a’’  sin“>,/  “  (Figure  A-4) 

2n  2 

and  g  (t)  is  as  before,  zero  for  |/j  >  t/2  (Figure  A-5),  the  convolution  integral  is 

pfz)  f  “  si"  ■  *>..  g  U)  r/i 

2n  ....  Hi  l) 


*H  IV  l  .iihi.  S mth/is.  </;///  ( 'nnwwnuatum.  John  Wilcv  ;nui  Son**,  p.  150. 
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or 


P(/)  =  f  sin  cu,/2(/  -  {,) 

2n  1  ' '  a),/2 (/  -  4) 


WM 


1 

-a 

‘X 

i 

<ux  <■> 

Figure  A-3.  Elementary  Gate  Function  Centered  On  Zero  Frequency 


This  bandwidth  limiting  causes  appreciable  distortion  in  the  original  pulse  shape  (Figure  A -6); 
however,  this  does  not  appear  to  degrade  the  resulting  pseudo-range  measurement  to  any  great 
degree.  The  GPS  receive!  cross-correlates  the  received  pulses  or  chips  with  a  locally  generated  code. 
The  result  of  the  cross-correlation  of  one  local  chip  with  the  received  bandlimited  chip  is  shown  in 
Figure  A-7.  This  correlation  curve  is  then  advanced  one-half  chip,  t/2,  and  retarded  one-half  chip. 
The  retarded  chip  is  then  subtracted  from  the  advanced  chip  defining  a  zero  crossing,  which  in¬ 
dicates  the  exact  time  of  maximum  correlation.  The  result  of  these  manipulations  is  shown  in  Figure 
A-8.  This  is  a  much  more  sensitive  method  of  obtaining  maximum  correlation  because  the  actual 
correlation  curve  (Figure  A-7)  has  a  rounded  top,  whereas  the  advanced  minus  retarded  curve  has  a 
well-defined  zero  crossing  due  to  the  steep  slope.  This  method  of  finding  maximum  correlation  is 
called  a  delay  lock  loop  * 


If  there  were  no  bandwidth  limiting,  the  chip  shown  in  Figure  A-5  would  be  received  with 
perfect  fidelity.  Consequently,  the  cross-correlation  function  would  have  a  triangular  shape  (Figure 
A-9)  instead  of  the  smooth  curve  of  Figure  A-7.  The  shape  of  the  curve  from  Figure  A-9 
corresponding  to  that  of  Figure  A-8  is  drawn  as  Figure  A- 10.  These  two  curves  are  more  closely  ap¬ 
proximated  if  the  bandwidth  of  the  receiver  is  increased.  This  result  is  illustrated  in  Figures  A-ll 
through  A-14  for  the  case  where  cu,  =  2rt/r.  and  in  Figures  A-15  through  A-I8  the  case  where  tu,  = 
3n/r. 


‘R  (  .  Dixon,  S/>rnnl  S/u’i  inim  Swtrim.  .Min  Wiles  anil  Sons,  p.  2<>S . 
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f  igure  A-4.  Frequency  Domain  Gale  Function  Transformed 
lo  the  Time  Domain,  te„  =  PI/TAU 
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F  igure  A-12.  Convolution  of  Momentary  Clare  With  Transformed  Clate.  ce,,  =  2  PI/TAU 
< Pulse  Shape  of  Received  Chips  l^iie  to  F  inite  Bandwidth) 
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The  single-frequency  pseudo-range  is  used  to  obtain  the  local  clock  correction.  As  was  il¬ 
lustrated  in  Figure  2  of  the  text,  this  correction,  tm  can  be  obtained  if  r,  (the  satellite  offset  from 
GPS)  and  t„  (the  satellite  offset  from  the  local  clock)  are  known.  The  predicted  values  for  t,  can  be 
obtained  by  decoding  the  satellite  navigation  message.*  The  corresponding  value  for  r0  can  be  found 
if  the  true  slant  range  to  the  satellite  is  known  at  the  time  of  each  observation.  A  reference  trajec¬ 
tory,  or  the  data  from  the  navigation  message,  can  be  used  to  determine  the  satellite  position.  At  the 
time  transmission,  let  this  be  r ,(i, ).  The  receiver  position  in  inertial  space  at  the  time  of  reception  is 
r  „(t2).  The  slant  range  r  is  the  difference 

r{l2)  =  r,(/,)  -  r,(i2) 

Therefore,  the  propagation  time  is  t,  =  r/c.  Propagation  delays  due  to  refraction  are  not  included 
in  this  calculation,  since  the  errors  they  contribute  are  less  than  the  1-jjs  accuracy  desired. 

The  pseudo-range  measurement  is  concluded  at  t2  when  the  satellite  epoch  is  received.  Conse¬ 
quently,  the  local  offset  will  be  time  tagged  /2,  which  equals  the  satellite  epoch  of  transmission  t, 
plus  the  correction  to  the  satellite  clock  t,  plus  the  propagation  delay  rr.  Therefore,  the  local  offset 
is 

t »Ui)  =  T„(r,  +  t,  +  rr )  =  Tp  -  Tr 
From  Figure  2,  the  correction  to  the  local  epoch  t,  must  be 

T,  =  T,  -  To  =  T,  -  Tp+  Tr 

where  t,  is  obtained  from  the  satellite  NAVDATA  message  and  is  calculated  using  the  station  and 
satellite  positions  at  the  times  of  reception  and  transmission,  respectively.  The  pseudo-range  obser¬ 
vation  is  Tp. 

The  data  t,  and  its  time  argument  t2  can  be  collected  and  fitted  by  a  polynomial  whose  coeffi¬ 
cients  were  obtained  using  least  squares.  The  polynomial  would  have  the  form 

t ,(/)  =  ^  aktk 

AsO 

Using  this,  the  local  time  is  corrected  by  putting  the  apparent  local  time  in  for  i  and  calculating  the 
correction  t,U).  Several  iterations  may  be  used  to  improve  the  accuracy. 


•Van  Oicremlonck.  Russell,  Kopit/ke.  Rirnhaum.  "The  CiPS  Navigation  Message",  Navigation,  Vol.  25,  No.  2.  Summer  1978. 
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This  simplified  development  establishes  a  worst  case  configuration  for  relative  satellite-earth 
station  motion.  It  is  useful  in  that  it  sets  limits  on  the  maximum  expected  Doppler  frequency  (veloc¬ 
ity)  and  Doppler  rate  (acceleration). 


The  worst  case  is  presumed  to  be  a  satellite  at  GPS  altitude  in  a  retrograde  equatorial  orbit 
(Figure  C-l).  No  GPS  satellite  is  ever  expected  to  be  in  such  an  orbit,  so  the  results  calculated  here 
will  determine  an  upper  limit  on  the  velocity  and  accelerations  along  the  slant  range  vector. 


Figure  C-l.  Retrograde  Equatorial  Orbit  Geometry 


From  the  law  of  cosines 
r  =  r,  +  r„-  2 rxrt  cos  Q 
Then 

2 r  dJL  =  2r,r„  sin  6  dl 
(it  dl 

Assume  that  oj,  and  a>,  are  constant. 

d®  =  CD  =  O),  +  CUf 

(it 
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(C-l) 


(C-2) 


(C-3) 


_  =  v  =  r>rKr~'  cd  sin  0 
dt 

^  =r<rK a,  fr  'cosfl  -  sin  0r'2  (C-4) 

dt  L  dt  dt  J 

Set  dv/dt  =  0  to  find  the  0.  which  gives  the  maximum  velocity: 


r 2  cos  6  =  f\/\.  sin 2  9 


This  can  be  rewritten  by  substituting  for  r1: 


r~  +  r‘  cos  0=1+  cos 1  9  (C-5) 

r.r. 

This  result  is  independent  of  ca.  Substituting  r,  =  26560  km,  rK  =  6378  km  results  in  a  value  for 
0  =  76°.  The  combined  angular  velocity 

fl)  =  2n  +  4n  _  6rr 
86164  86164  86164 


Evaluating  Equation  C-3  at  9  -  76°  gives  the  maximum  velocity 


r'r*u  Sln  _ =  1395  m/s  (C-6) 

[r,J  +  rl  -  Irj,.  cos  0] ,/J 


At  the  two  frequencies,  L,  =  1575.42  MHz  and  L,  =  1227.6  MHz,  the  Doppler  shift  is  approx¬ 
imately  7326  Hz  on  L,  and  5708  Hz  on  l.t.  The  acceleration  is  obtained  by  substituting  Equation 
C-3  into  Equation  C-4: 


dv 

dt 

dv 

dt 


=  r,rBa)[r"'co  cos  0  -  sin2  0  /•'V,rtcu) 
=  r,rsaj2r  '[cos  0  -  r’V.r,  sin  2  0] 


By  inspection,  it  is  evident  that  dv/dt  will  be  greatest  when  0  =  0.  The  second  term  in  brackets  is  zero 
when  0  =  Oand  is  greaterthan  zero  for  any  other  0.  Also,  Figure  C-l  showsthat  rwill  beaminimum 
at  0  =  0,  which  makes  r  '  and,  consequently,  dv/dt  a  maximum.  At  0  =  0,  Equation  C-4  becomes 


dv 

Jt 


r,r„ ou 


|/\2  +  r‘  -  2 r.r,]ul 


=  0.402  m/s2 
56 


t1!'. 


The  maximum  Doppler  rate  is  also  at  the  angle  of  maximum  acceleration.  At  the  two  GPS  frequen¬ 
cies,  these  are  2.1  Hz/s  on  L{  and  1.6  Hz/s  on  L2. 


The  error  in  the  Doppler  count  over  1  min  due  to  the  use  of  the  average  frequency  can  be 
estimated  using  these  rates.  With  an  L ,  rate  of  2.1  Hz/s,  the  average  frequency  can  be  in  error  by 
about  63  Hz  at  each  end  of  the  60-s  count  interval.  Assuming  a  value  for  the  fractional  interval  equal 
to  a  full  cycle  (using  the  zero  Doppler  frequency  28750  Hz),  the  count  will  be  in  error  by  63/28,750 
=  0.002.  This  is  equivalent  to  a  0.04-cm  error  in  the  delta  range  measurement  at  L, .  The 
corresponding  error  at  L2  is  0.03  cm. 
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